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We discuss a protocol for growing states with topological order in interacting many-body systems 
using a sequence of flux quanta and particle insertion. We first consider a simple toy model, the 
superlattice Bose Hubbard model, to explain all required ingredients. Our protocol is then applied 
to fractional quantum Hall systems in both, continuum and lattice. We investigate in particular how 
the fidelity, with which a topologically ordered state can be grown, scales with increasing particle 
number N. For small systems exact diagonalization methods are used. To treat large systems with 
many particles, we introduce an effective model based on the composite fermion description of the 
fractional quantum Hall effect. This model also allows to take into account the effects of dispersive 
bands and edges in the system, which will be discussed in detail. 


I. INTRODUCTION 

Since the discovery of the integer quantum Hall ef¬ 
fect (IQHE) in 1980 [I| and two years later the frac¬ 
tional quantum Hall effect (FQHE) 0, the interest in 
states exhibiting topological order has increased tremen¬ 
dously. Due to their robustness against local disorder, 
these states are for instance well suited for metrology 
[l| . Another interesting aspect of such interacting many- 
body systems is that they host states with anyonic ex¬ 
citations il. In this context, non-Abelian anyons are 
particularly exciting because they can be employed to 
build a topological quantum computer Q , if they can be 
coherently manipulated. 

For many years, solid state systems with electrons ap¬ 
peared to be the only candidates to realize exotic many- 
body states, e.g. in the FQHE. However, their small in¬ 
trinsic length scales make coherent control difficult. On 
the other hand, the rapid development of ultracold gases 
and photonic systems is a promising route to implement 
various Hamiltonians known from the solid state com¬ 
munity, with full coherent control. With current state 
of the art technology, noninteracting tOTological states 
have been observed in ultracold gases [7Hl0l| as well as 
photonic systems [ll|, [ 13 . 

Unlike in the solid state context, efficient cooling mech¬ 
anisms below the many-body gap are still lacking in sys¬ 
tems of interacting atoms and photons, which makes the 
preparation of topological states challenging. In this pa¬ 
per we develop an alternative dynamical scheme for the 
efficient preparation of highly correlated ground states 
with topological order. It can be implemented with ultra¬ 
cold atoms or photons and exploits the coherent control 
techniques available in these systems. 

To solve the cooling problem, different approaches were 
discussed previously. In [T^ it was suggested to pump 
the sites of a dissipative lattice system coherently. This 
scheme relies on a direct A^-photon transition and thus 
works only for small systems. Moreover, the state pre¬ 
pared is in a superposition of different particle numbers. 
Another proposal suggested to rapidly refill hole exci¬ 
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FIG. 1. (Color online) Simulation of the effective composite 
fermion (CF) lattice model. The figure shows the CF den¬ 
sity p* of a lattice with 11 x 11 sites and effective fiux per 
plaquette a* = 0.1. The basic ingredients of the protocol are 
(i) fiux insertion, which creates a CF quasi-hole excitation 
and (ii) a coherent CF pump, which refills the hole excita¬ 
tion. Additionally, we include a quasi-hole trapping potential 
(ph = 1'^*) and absorbing boundaries (yndge = 0.0025J*) to 
account for the effects of dispersive bands and edge states. 


tations due to local loss in a lattice system. This protocol 
stabilizes the ground state dynamically. 

In reference [l^ . we proposed a scheme to grow the 
highly correlated Laughlin states. Here, we discuss the 
general ingredients of the protocol, which may also be 
used to grow other topological ordered states in interact¬ 
ing many-body systems. To do so, we consider a simple 
toy model explaining all necessary ingredients. The main 
ideas of the protocol can be summarized as follows. In the 
first step, a topologically protected and quantized Thou- 
less pump M is used to create a localized hole excitation 
in the system. In the second step the hole is refilled using 
a coherent pump. Local repulsive interactions between the 
particles are employed to provide a blockade mechanism 
such that only a single particle is inserted. In the case of 
the FQHE, we analyze the performance of the protocol 
in detail and explain how the fidelity scales with the par¬ 
ticle number N. In order to treat large fractional Chern 
insulator systems with many particles, we introduce an 
effective model based on the composite fermion (CF) de¬ 
scription of the FQHE. Within this model, we are able 
to demonstrate, that the growing scheme still works de¬ 
spite the presence of gapless edge states and dispersive 
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bands in large lattice systems. As shown in Fig. [T] we 
reach a homogeneous density in the bulk with an average 
CF magnetic filling factor ly* ~ 0.9, close to the optimal 
value ly* = 1. 

The paper is structured as follows. We start by de¬ 
scribing the growing scheme using a simple toy model 
in section im In section uni we discuss the protocol in 
a FQH system. Moreover, we discuss the performance 
of the protocol. The last section IIVI considers fractional 
Chern insulators. There, we discuss edge and band dis¬ 
persion effects based on our effective CF model. 


II. TOPOLOGICAL GROWING SCHEME 

In this section, we discuss the basic idea how states 
with topological order can be grown for the one dimen¬ 
sional superlattice Bose Hubbard model (SLBHM) as a 
toy model. This sets the stage for the following discus¬ 
sions of two dimensional systems with topological order. 


A. Toy model 

The SLBHM can be used to implement a quantized 
pump which can be related to a nontrivial topo¬ 

logical invariant. The Hamiltonian of the SLBHM is 

a = - Ji ^ a^-hj + h.c. - -^2 ^ + h.c. 

3 3 

3 3 

where aj{bj) annihilates a boson on lattice site A {B) in 
the jth unit cell (see Fig. [2^). Note, that we consider 
here a semi-infinite system with open boundary on the 
left. The hopping elements are denoted by Ji, J 2 , while 
S is an onsite potential shift acting on lattice sites B. 
Furthermore, we include onsite interactions U for more 
than one particle per site. 

In the case of vanishing S and disregarding interactions, 
the different hopping amplitudes Ji ^ J 2 determine the 
dimerization of the two sites A and B in the unit cell. 
This leads to a two-band model with a band gap Agap 
and bandwidth AE determined by the ratio J 2 /Ji- In 
the limit J2/J1 0 the flatness ratio Egap/AE tends 

to infinity and the two bands become nondispersive. At 
filling p = 1/2 and with interactions, the ground state 
of the system is a Mott insulator (MI) with many-body 
gap A or a superfluid (SF) depending on the specific 
parameters J2/J1 and U (isl - l^ . Specifically, for large 
interactions U the model can be mapped to free fermions. 
The resulting ground state at p = 1/2 is incompressible 
whenever Ji ^ J 2 or |(5| 0. For finite interaction U, 

the SF region is extended in parameter space (Ji — J 2 , <5) 
as depicted in Fig. [^Jo. 



FIG. 2. (Color online) (a) Semi-infinite chain of coupled 
dimers, consisting of two sites A, B. Initially, the MI ground 
state is prepared with filling p = 1/2. (b) Schematic phase 
diagram of the SLBHM in the regime of intermediate inter¬ 
actions U. Following adiabatically the loop 1-5 corresponds 
to a Thouless pump: Encircling the SF region, all particles in 
the MI are shifted one dimer to the right. 

B. Protocol 

In the following, we illustrate the main steps of our 
growing scheme by showing how the MI ground state of 
this system can be grown. We note that there are more 
practical ways to prepare the ground state, of course, but 
the present protocol can be generalized to more complex 
systems, including states in the FQHE presented in sec¬ 
tion mil and IIVI 

We will start to discuss the protocol in the nondisper¬ 
sive limit J2/J1 —^ 0, where hole excitations will remain 
localized. The effect of dispersive bands in the case of 
finite J1/J2 will be considered later. Let us assume, that 
the N particle MI ground state of the SLBHM with fill¬ 
ing p = 1/2 is already prepared in a finite region of the 
lattice. We now show how the state with A -|-1 particles 
can be grown. 

1. Topological pump 

First, a topological Thouless pump [I^ is used to cre¬ 
ate a hole excitation localized at the edge of the system. 
This process is related to Laughlin’s argument of flux in¬ 
sertion [23 in the quantum Hall effect (QHE): Inserting 
one flux quantum c/o leads to a quantized particle trans¬ 
port, which is the origin of the quantized Hall conduc¬ 
tance axy In the SLBHM, the process of flux insertion 
is shown in Fig. [2 )d. By adiabatically changing the pa¬ 
rameters Ji — J 2 and 5 in time we encircle the critical 
SF region. After a full cycle of period Tfj, all particles are 
shifted one dimer to the right, leaving a hole excitation on 
the left side of the system. The possibility to construct 
such a cyclic pump is deeply related to the underlying 
topological order of the model. 

2. Coherent pump 

In the following step, we replenish the hole excitation 
by a single boson. To this end, we coherently couple a 
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FIG. 3. (Color online) (a) Coherent pump couples a reser¬ 
voir of bosons to the left-most dimer, (b) Due to a blockade 
mechanism, only a single particle is pumped into the system, 
fleff is the many-body Rabi-frequency, reduced by a Franck- 
Condon factor and A^o the number of bosons in the left-most 
dimer. 


reservoir of particles to the first lattice site of the system 
as shown in Fig. [3^. The coherent pump can be described 
by the Hamiltonian 

Tin =F!aJe-“‘-fh.c., (1) 

where 17 is the Rabi- and uj the driving frequency. The 
driving frequency uj is chosen to be resonant with the 
lowest band. Hence, if |17| is much smaller than the 
single-particle band gap F^gap, particles can only be added 
in the lowest band. The corresponding Rabi-frequency 
|Heff| < 1171 is then however reduced by a Franck-Condon 
factor. 

In general, the coherent pump ([I]) creates a coherent 
superposition state of Nq bosons in the left-most Wannier 
orbital of the lowest band. To ensure that, at most, one 
boson is added, we employ a blockade mechanism, see 
Fig. [3 ]d. For sufficiently strong interactions and weak 
enough driving, 


C/»|17eff| (2) 

the transition to the Nq = 2 particle state is detuned by 
the interaction energy U. 

For |17eff| <C A, at maximum one boson can be added 
into the lowest Bloch band. To replenish the hole exci¬ 
tation by precisely one boson at a time, one can use a 
TT-pulse of duration Tq = 7r/217eff. Although the 7 r-pulse 
is less robust to errors than e.g. an adiabatic sweep, we 
choose it because of its speed. This, we believe, is crucial 
to overcome linear losses present in realistic systems. 

Starting from complete vacuum, the sequence de¬ 
scribed above can be repeated to grow the MI ground 
state with N particles. Next, we discuss extensions of 
the protocol, required when the bands are dispersive or 
the system is finite. 



FIG. 4. (Color online) Due to the process of flux insertion 
in a finite system, the higher band of the SLBHM will be 
occupied. While a state at the right edge of the system will 
be transferred into the upper band, a state from the left upper 
band, will be transferred to the lower band. 


bandwidth AE of the lowest band sets the timescale for 
the dispersion. Firstly, the localized quasi-hole excita¬ 
tion created using the topological pump, disperses into 
the bulk of the system. Since the coherent pump cou¬ 
ples locally to the left-most dimer, the hole cannot be 
refilled efficiently with a single boson. More importantly. 


the number fluctuations AN{t) 


{N^)-{Ny /{N) 


will increase. To prevent the dynamics of the quasi-hole 
excitation, a quasi-hole trapping potential 



(3) 


acting on the left-most dimer can be used. The strength 
of the trap gh should be weak enough not to destroy the 
topological order, but strong enough to trap the quasi¬ 
hole excitation. Concretely we require 


AA < 5h "C a. 


(4) 


Secondly, the diffusion of bulk particles from the right 
edge of the system into vacuum lets the MI melt. Vice- 
versa, the diffusion of holes from the vacuum into the 
bulk at the right edge of the system makes the state 
compressible. For fast pump rates of the protocol, this 
effect can be neglected. However, if the bandwidth AE 
is the dominant contribution, we suggest to use a har¬ 
monic trapping potential. The potential should be weak 
enough to avoid localization effects (see [l^). Therefore, 
in local density approximation, the bulk of the system 
will remain incompressible. 

In finite systems with open boundaries, the flux in¬ 
sertion, as illustrated in Fig. IH connects states from 
the lower and the upper band at the edges. As long as 
the upper band is empty, the topological pump creates 
a hole excitation in the lower band. However, once the 
right boundary is reached, particles will be excited to the 
upper band. To avoid the high energy excitations in the 
protocol, sufficiently strong losses 7 Edge localized at the 
right boundary of the system can be introduced. 


III. FRACTIONAL QUANTUM HALL STATES 


3. Dispersive bands & finite systems 

Dispersive bands J 1 /J 2 7 ^ 0 lead to intrinsic dynam¬ 
ics of the particles and thus of the hole excitations. The 


We now apply the growing scheme to fractional quan¬ 
tum Hall (FQH) systems. In this case, the bands - i.e. 
Landau levels (LLs) - are nondispersive and we assume 
an infinite system. 
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A. Model 

We consider a FQH model of bosons in the two dimen¬ 
sional disc geometry. The magnetic field can be iirrple- 
mented using e.g. artificial gauge fields. Moreover, we 
consider a contact interaction between the particles with 
strength g^. In second quantized form, the Hamiltonian 
reads 

= y d^r [p - Af 

+ ^90 j d^r (5) 

where creates a boson at the position r. The first 

term in eqn. ([5]) includes the magnetic field in miirimal 
coupling using the vector potential A = B/2{—y,x,0). 
In this symmetric gauge, the total angular momentum 
Lz is a conserved quantity, [H, Lz] = 0. 

For later purposes we express the held operators 
in terms of the bosonic operators which create a par¬ 
ticle in the orbital of the nth LL with angular momentum 
£. Here, n and i are integers with n, £-|-n > 0. We obtain 

i’ir) =^Vn,e{r)bn,e, (6) 

wher e ij„j( r) are the single particle wavefunctions, see 

e.g. dlil. 

B. Laughlin States and Excitations - CF Pictnre 

It was shown in [26l - [30l| . that the Laughlin (LN) wave- 
function at magnetic Filing z/ = 1/2 is the exact zero en¬ 
ergy ground state of the bosonic FQH Hamiltonian ([5|). 
The Filing factor z/ = N/N^ is deFned as the ratio be¬ 
tween particle number N and number of Fux quanta fV^. 
In terms of the complex coordinates Zj = Xj — iyj of the 
jth particle, the LN wavefunction is 

|LN,7V) = 4-ln = A4.N (7) 

j<k 

where A/ln is a normalization constant. The magnetic 
length £b = a/A/B depends only on the magnetic Feld 
B. 

The zero energy excitations of the LN wavefunction 
are quasi-holes described by the wavefunction 4/^^, = 
A/qh rij •Zj'I'ln- tti quasi-holes, located in the center, are 
described by the m quasi-hole wavefunction 

Imqh, A) = d/mqh — A/rnqh (n^j) (8) 

3 

The LN state |LN, A^) and the quasi-hole state |mqh, A^) 
carry different total angular momentum Lz, 

L,(|LN,iV)) =iV(iV-l) (9) 

Lz{\m(^ii,N))=N{N-l) + mN. (10) 


Both wavefunctions can be understood in the more 
general framework of composite fermions (CFs) miii- 
1841 . For the z^ = 1/2 case, the LN wavefunction ([71) can 
be decomposed into 

'I'ln = A/ln W _ i .^3 ~ (H) 

j<k 

Besides the Jastrow factor ~ attaching 

one Fux quantum to each boson (see Fig. [5 )d), a 
fermionic wavefunction appears. The Jastrow 

factor screens the interactions between the particles and 
leads to a reduced magnetic Feld B* = B/2 seen by 
the CFs. As in the IQHE, their wavefunction is given 
by a Slater determinant of single particle orbitals Fil¬ 
ing the lowest CF-LL v* = 1, i.e. = Y\j<ki ^3 ~ 

Zk) exp |ZjP/4£|^. The FQHE can be interpreted 

as an IQHE of noninteracting CEs in a reduced magnetic 
Feld B*. The CE picture in the context of the EQHE is 
powerful in describing all fractions in the Jain sequence 
and moreover describes the quasi-hole excitations 
and their counting correctly. However, the CF theory 
does not reproduce the correct Laughlin gap Aln- Here, 
a microscopic theory of the full interacting many-body 
problem is necessary. Also in order to explain other than 
the main sequence of FQH states, interactions between 
CFs need to be taken into account. 


C. Protocol 

1. Topological pump - flux insertion 

In the toy model III Al we used a topological Thouless 
pump to shift the particles to the right. In the context of 
quantum Hall physics, this corresponds to Laughlin’s ar¬ 
gument of flux insertion [^ . which was used to explain 
the quantization of the Hall conductivity. The idea is 
to insert locally in the center of the system Fux quanta 
(j)Q in time which produces an outwards Hall current 
jr ~ Ch dt4> in radial direction. The quantization of the 
Hall current is related to the nontrivial Chern number 
Ch of the system. The process of inserting flux quanta 
increases the total angular momentum of the state. A 
more detailed discussion of the process in the noninter¬ 
acting case can be found in appendix lAl 

Starting from the z/ = 1/2 LN state |LN, A) with N 
particles, we insert two flux quanta 2 x c/o in time in 
the center of the system. In this way, we create a two- 
quasi-hole excitation |2qh, A). The two-quasi-hole state 
|2qh, A) already lies in the same angular momentum Lz 
sector as the LN state |LN, N + 1) with A J- 1 particles 
(compare eqns. dHl), (flUl) '). 
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FIG. 5. (Color online) (a) Growing scheme for the v = 1/2 LN state in the continuum. The many-body eigenstates of the FQH 
Hamiltonian between which transitions are driven, are shown as a function of the total conserved angular momentum Lz- 
Starting from the vacuum |0), the LN state with N = 2 bosons is grown. After adding one boson in the central orbital, two 
flux quanta are adiabatically inserted and the resulting hole excitation is refilled by a second boson. This can be interpreted 
as replenishment of the system by a composite fermion, hence keeping intact the topological order of the system. The lowest 
line shows the schematics of the prepared states after every step. (Horizontal lines denote orbitals; This picture is valid in the 
thin-torus limit or after using the patterns as partitions determining certain Jack polynomials [^.) (b) Flux attachment: A 
single boson together with one flux quantum form a GF. (c),(d) Level spectrum in terms of CFs. (c) Goherent pump process 
with Rabi-frequency H. The lowest GF-LL is filled, corresponding to GF filling factor v* = 1. (d) Flux insertion process 
increases angular momentum I by one. 


2. Coherent pump 

In the next step, we refill the quasi-hole excitation us¬ 
ing a coherent pump, which couples a reservoir of bosons 
to the hole excitation. As this is done in the center, no 
angular momentum is transferred. To this end, we sup¬ 
plement the FQH Hamiltonian ([5]) by a term 

■^0 = y d^r g(f)e““*(5^(r)'0'l'(r)-F h.c. (12) 

While g{t) denotes the strength of the coherent pump, 
the driving frequency oj = is chosen resonant with 
the lowest LL (LLL). In the regime of large magnetic 
fields and low temperatures, it is sufficient to project eqn. 
(fT^ to the LLL resulting in 

'Plll^p’Plll = g -F h.c., (13) 

where we defined the Rabi-frequency O = gj •\/27r.^g. 
This corresponds to Ileff in the toy model. 

Since the coherent pump does not change the angu¬ 
lar momentum sector, we obtain a similar block¬ 
ade mechanism as in the toy model (see Fig. with 
Lz. = 0, 2). In the corresponding sector, there is only one 
zero-interaction energy eigenstate, the N + 1 LN state. 
The energy offset to any other states in the {N -F 1)- 
particle manifold is of the order of Haldanes zeroth pseu¬ 
dopotential Vq [ 3 ^. We require H <C WcVb to avoid 
excitations to high energy states. 

To insert a single particle, we use a 7r-pulse of dura¬ 
tion Tq = Tr/2Q^^\ Here, the bare Rabi-frequency fl is 
reduced by a Franck-Condon factor 

HW/L!=(LN,iV|^_o|2qh,iV-l), (14) 



1 /N 


FIG. 6. (Color online) Franck-Condon factors as de¬ 

fined in eqn. d, for different particle numbers N calculated 
using exact diagonalization (ED) in a system with contact in¬ 
teraction for up to A = 7 particles. A linear fit is used to 
extrapolate A —>■ 00 . 


which accounts for the overlap between the initial and 
final state. In Fig. [6]we show the Franck-Condon factors 
for different particle number N and extrapolate to A —F 
00 . To this end, we calculated the LN wavefunction using 
exact diagonalization [s^ . A linear fit suggests that even 
in the case A —>■ 00 the Franck-Condon factor does not 
vanish with /H ~ 0.7. 


3. CF picture 

The first few steps of the protocol are illustrated in 
Fig. [S^. Another description of the protocol, which pro¬ 
vides a much simpler physical picture uses the CF pic- 
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ture. Firstly, due to local flux insertion (see Fig. [5jl), we 
generate a quasi-hole excitation. Secondly, the hole exci¬ 
tation is refilled by an effective coherent CF pump, for a 
particle together with one flux quantum. Analogously to 
the blockade mechanism, one CF refills the empty £ = 0 
orbital (see Fig. [3:). By subsequent repetition, we grow 
the filling i/* = 1 integer quantum Hall state of noninter¬ 
acting CFs in a reduced magnetic field H*. 

The CF picture provides an explanation, why it is pos¬ 
sible to grow highly correlated LN states. Since the CF 
wavefunction is a separable Slater determinant 

of noninteracting CFs, we can grow the highly correlated 
LN state by adding CFs one by one to the system. Note 
that, although the wavefunction is a separable product 
state of CFs occupying Wannier functions (or Landau 
level orbitals), it has non-local topological order. This 
is because Wannier functions (or orbitals) of a single LL 
are non-local [s^- On first glance, this may look con¬ 
tradicting as the dynamical growing scheme presented in 
this paper involves only local processes. However, the 
non-local topological order in the system is generated by 
the Thouless pump (flux insertion). 


D. Performance 


For interesting applications such as measuring the 
braiding statistics of elementary excitations, the LN state 
has to be prepared with high accuracy. To quantify the 
quality of the scheme, we calculate the fidelity of be¬ 
ing in the LN state with N particles after N steps of the 
protocol. The fidelity is defined as 

J-^ = |(«'jv|LN,Ar)|, (15) 

where I'f'Ar) is the state after N steps of the protocol. 
In the calculation of the fidelity we include besides 
nonadiabatic transitions in the flux insertion and coher¬ 
ent pump process, a particle loss rate 7 . This is impor¬ 
tant, since decay usually plays a crucial role, in particular 
in photonic systems. 

The total time for a full cycle (flux insertion time 
coherent pump time Tq) is T = + Tq. As will be 

discussed in detail later, we find for the fidelity pertur- 
batively in the limit 7 T, (Aln?!^)”^, (AlnIo)”^ <C 1 

.r„.exp(-iv(i-,r(« + i) + ^^)). , 16 ) 

We see, that the losses and the nonadiabatic transitions 
contribute competitively in eqn. (HU). While it is favor¬ 
able to run the protocol as fast as possible to avoid losses, 
the adiabaticity requires long time scales T. 

For fixed fidelity = 1 — e, e<C 1, we calculate the 
maximal number of particles A^max which can be grown 
with optimal period Topt. To leading order in JV, we 


approximately obtain 

/A \ 2/5 

1.365£3/5 j (17) 

Topt 1.431 (£7)"'^' ' ■ (18) 

The LN state with fildelity J-jy = 1 — £ can be grown in 
time To = A^maxTopt ~ Arm'ix£“^/^AAr/ALN scaling only 
slightly faster than linear with particle number A^max- 
Note, that the fidelity decreases with increasing time 
t > Tq. In the following, we discuss the different con¬ 
tributions to the fidelity (|T 6 )l in detail. 


1. Loss 

We assume a particle loss rate 7 1/T, i.e. after one 

full period T the probability of loosing a particle is small. 
The probability of a single decay process after N steps 
of the protocol is then given by 

Ty = 1 - exp ^-yT^n^ ~ i 7 TA(A-h 1). (19) 

To leading order in N, the probability of loosing one 
particle increases quadratically with the particle number 

N. 


2. Flux insertion 


The process of flux insertion, as described in IHI Cl in¬ 
troduces one flux quantum c/jq in the center of the sys¬ 
tem. For simplicity, we assume the flux = (po t/T^ 
to change linear in time t. In the fully adiabatic proto¬ 
col, the angular momentum £ is therefore increased by 
one without coupling different LLs n,n'. We calculate in 
the noninteracting case the nonadiabatic coupling matrix 
element k = {n',£\ — id(,\n,£) between different LLs. To 
estimate the scaling of the probability for excitation 
of higher LLs, we consider the coupling between the LLs 
n and n-l-1. In the regime of long times T^, we determine 
perturbatively (see appendix iBl) . Using the Laughlin 
gap Aln, we obtain 


“ (AlnT^)2 ’ 


( 20 ) 


where k is a nonuniversal coupling constant, 
tantly, nonadiabatic transitions to higher LLs 
1/(AlnT,^)^. 


Impor- 
scale as 


3. Coherent pump 

The coherent pump (fT3) . as discussed in IHI Cl couples 
in zeroth order in H/Aln only the two-quasi-hole state 










7 


1.5 

1 

CTat 

0.5 

0 


- ^0- - ^ ' 


0 ED 

- - quadr. fit 


0 0.2 


0.4 0.6 0.8 

\/N 


1 


FIG. 7. (Color online) Scaling of the nonuniversal factor crjv 
from eqn. mi with particle number N. We interpolate N 
oo with a quadratic fit. 


effects in our investigation. Moreover, we consider finite 
systems now, where edge states are present. 


A. Hofstadter Hubbard Model 

We consider a two dimensional lattice described in a 
tight-binding model with nearest neighbor hopping terms 
J. The Peierls phases are chosen to mimic a magnetic 
field in Landau gauge. Additionally, we consider on¬ 
site interactions U resulting in the Hofstadter Hubbard 
Hamiltonian 

a = ^ + h.C. 

x,y 

A U/‘^ .y^t:,y<lx,y (23) 

x,y 


|2qh, N) to the Laughlin state |LN, -|- 1). In first order 
perturbation theory, couplings to excited states in the 
N — 1, N, N + 1, N + 2 particle sectors are relevant. We 
obtain for the probability of nonadiabatic excitations 


p ^ 

(AlnTo)2’ 


( 21 ) 


where we used the 7r-pulse time Tq = 7r/2H. The nonuni¬ 
versal factor (7AT includes Franck-Condon factors and ex¬ 
citation energies of undesired states weighted by 

the Laughlin gap Aln (see appendix [Cl. In Fig. (TIctat 
is calculated for various particle numbers N using ex¬ 
act diagonalization. For A —> oo we extrapolate with a 
quadratic fit — 1-4. 

To conclude, we obtain the estimate for the fidelity in 
eqn. dm, summarizing the contributions of flux insertion 
and coherent pump Pq for fixed ratio Tq/T^. Assum¬ 
ing, the coefficients k, ctat to depend only slightly on the 
particle number N, we define the factor Aat containing 
all constant contributions. Then, using 

Pn = exp (22) 


we obtain eqn. dm- 


IV. LATTICE AND FRACTIONAL CHERN 
INSULATORS 


The X, y coordinates are measured in units of the lat¬ 
tice constant a. We use the operators to denote 
the bosonic annihilation of a particle at site (x, y). The 
flux per plaquette a (in units of the flux quantum) ac¬ 
counts for the magnetic field penetrating the two dimen¬ 
sional lattice. In particular, it sets the magnetic length 
= a/^/2JTa, which describes the extend of the cy¬ 
clotron orbits. 

In recent experiments, the realization of the Hofstadter 
model was shown in a photonic system 11 ll a s well as ul¬ 
tracold gases [3, Hi- The photonic systein[ll| implements 
a tight-binding model using ringresonators in the opti¬ 
cal wavelength regime. The tunneling is achieved using 
off-resonant ringresonators, which have a different length 
for hopping forward and backward and the Peierls phases 
are determined by the optical path difference. In exper¬ 
iments with ultracold gases [3, Q, standing waves are 
used to create a two dimensional optical lattice. Due 
to a strong field gradient along one direction, the par¬ 
ticles are localized. Using laser-assisted tunneling tech¬ 
niques, the hopping elements are implemented with an 
additional Peierls phase. In the experiments, the Peierls 
phase can be engineered to realize e.g. a flux per plaque¬ 
tte of a = 1/4. 

Note, that the continuum limit a —^ 0, where the mag¬ 
netic length is much larger than the lattice constant 
a, corresponds to the FQH model dSI. 


B. Ground State and Excitations 


Lattice systems are promising candidates to realize 
FQH physics, since the magnetic helds realized in these 
systems are very strong, such that low magnetic filling 
factors v < 1 can be achieved while keeping a sufficiently 
large density of particles. However, in this case the mag¬ 
netic length is and the lattice constant a become of com¬ 
parable size. Therefore lattice effects, such as dispersive 
bands, play a crucial role. Now, we introduce an effec¬ 
tive CF lattice model, which allows us to include these 


We summarize the properties of the ground state of the 
Hofstadter Hubbard model (1^ for bosons with magnetic 
filling factor i/ = 1/2 in two different geometries. In the 
case of a torus discussed in [H, [s^, the exact ground 
state was compared to the Laughlin state 0 projected 
on a lattice for different flux per plaquette a. It was 
found, that the Laughlin wavefunction provides a good 
description up to a ~ 0.2. Moreover, the many-body 
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FIG. 8. (Color online) a) Flux insertion on a Ceo buckyball 
geometry, b) Many-body spectrum i?(0) on a Ceo buckyball 
starting with N — 3 particles and = 4 flux quanta and 
introducing two more flux quanta 2 x 0o. 


Chern number remains Ch = 1/2 until the flux per pla- 
quette reaches a critical value of a ~ 0.4. The size of the 
Laughlin gap Aln depends on the parameters a and U. 
For instance, at a ~ 0.1, the Laughlin gap saturates at 
Aln — O.IJ for large U ^ J. 

In [l^, we numerically analyzed the ground state 
in a spherical geometry using a buckyball lattice with 
Ns = 60 sites (see Fig. [5^). To realize a filling z/ = 1/2 
in the continuum on a sphere, the relation between par¬ 
ticle number N and flux quanta iV^ is iV^ = 2{N — 1). 
We identify the topological order of the ground state by 
inserting two flux quanta <00 in a system with N = 3 
particles and = 4 flux quanta. The many-body spec¬ 
trum during flux insertion is shown in Fig. [5 Jd. We 
observe a single incompressible ground state with many- 
body gap A ~ O.IJ for U = lOJ, similar to the case 
on a torus. Importantly, we obtain the correct count¬ 
ing of nearly degenerate quasi-holes states after inserting 
one and two flux quanta as expected from the continuum 
limit on a sphere. In both lattice cases we expect the 
ground state to be in the same topological universality 
class as the LN state in the continuum. 


C. Effective CF Lattice Model 

Before we introduce the effective CF lattice model, we 
briefly summarize the results of the exact simulations in 
reference [l^. We implemented the protocol on a Ceo 
lattice up to A = 3 particles. Additionally, we used a 
similar quasi-hole trapping potential as in eqn. ([3|). The 
numerical results show, that the state prepared af¬ 

ter three steps of the protocol is close to the LN ground 
state |gs) with N = 3 particles. Moreover, as expected 
from the blockade mechanism in the coherent pump, the 
number fluctuations AN(t) are small throughout the pro¬ 
tocol. 

To discuss edge effects, larger systems with many par¬ 
ticles are required, such that bulk and edge states can 
be distinguished. However, using exact diagonalization 
only a few particles are feasible on large lattice systems 
with Ns > 60 sites. To simulate such a system, we need 


an effective model describing the low energy dynamics of 
the Hamiltonian 

In the spirit of the composite fermion principle, we as¬ 
sume a model of noninteracting composite fermions on 
a lattice. By attaching one flux quantum to a boson, 
as shown in Fig. [SJd, an (approximately) noninteract¬ 
ing composite fermion is formed in a reduced magnetic 
field. The simplest tight-binding model with a reduced 
magnetic field, a* = a/2, only considers nearest neigh¬ 
bor hopping elements J*. Under these assumptions, the 
effective model is 

HcF = — J \^x+l,y^x,y^ ^ + C^x,y+Nx,y + h.C. , 

x,y 

(24) 

where cl, y creates a composite fermion at site {x, y). The 
only free parameter J* in this model determines the time 
scale of the dynamics. By comparing the bandwidth of 
the full many-body spectrum of interacting bosons to the 
free CF single-particle spectrum on a Ceo buckyball at 
N^ = 6 flux quanta, we find J* ~ J. 

Now, we conjecture that the effective CF lattice model 
describes correctly the low energy dynamics. Yet, there 
is no proof that the CF theory is applicable in a lattice. 
However, there are several hints that the essential physics 
can still be understood in terms of CFs. First of all, the 
CF picture in a lattice captures the correct counting of 
quasi-hole excitations [d^. Thus, the low energy exci¬ 
tations should be described correctly. Moreover, the LN 
wavefunction is a special case of the more general CF the¬ 
ory, when the lowest CF-LL n* = 0 is filled. As shown in 
the LN wavefunction projected on a lattice describes 
the many-body ground state on a lattice accurately up 
to relatively high flux per plaquette a ~ 0.2. Therefore, 
we limit the flux per plaquette in our effective model to 
a* < 0.1. Finally, CF states from the Jain sequence [sj, 
other than the LN state, have been identified in a lattice 
for small flux per plaquette. Therefore, we expect 
that this effective model captures the essential physics of 
the original many-body model in the lowest Chern band 
(LChB), including the dynamics of the hole excitations. 
As this model is supposed to describe only the low energy 
regime, excitations to higher Chern bands (HChBs) are 
not expected to be described correctly. 


D. Numerical Results 

1. Numerical implementation of the protocol 

The protocol is implemented in both, the spherical 
geometry on a Ceo buckyball as in [l^ as well as the 
square lattice with open boundary conditions. We use 
the method explained in [l5l| to insert flux quanta locally. 

Unlike the coherent boson pump dn, we model the 
coherent CF pump by coupling a CF reservoir mode to 
the central site of the system. The reservoir mode is 
refilled in each step of the protocol. Therefore, we insert 
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at most one CF per cycle. This is similar to the blockade 
mechanism for bosons and therefore only justified in the 
limit of large interactions U ^ J, where corrections scale 
as (O/Aln)^ (see eqn. (HD)). 

Moreover, we implement the quasi-hole trapping po¬ 
tential by including an onsite potential gh on the central 
site. 

As explained for the toy model, we introduce loss chan¬ 
nels at the boundaries of the system to prevent high en¬ 
ergy excitations. Then, we calculate the time evolution 
of the correlation matrix elements {c\.^yCxi^y')■ by solv¬ 
ing the corresponding master equation in Lindblad form 
{h=\) 

1 ' . , 

dtP = P] + 2 ^ / ‘^^x,yPlx^y ~ {^x,y^x,yj P} (25) 

x,y 

numerically. Absorbing boundary conditions are de¬ 
scribed by the jump operators lx,y = ^/TidgeCx.y with 
loss rate 7 Edge and we restrict the sum to the edge of 
the system. 

Note, that the loss of a CF is a loss of both, a flux 
quantum and the boson it was attached to. Therefore, we 
only allow this loss term at the edge of the system, where 
the meaning of a free flux quantum without a particle is 
obsolete. 


2. Performance 

To investigate the performance of the protocol, we use 
the spherical geometry without edges. The bandwidth 
AE/J of the Ceo buckyball with up to = 6 flux 
quanta is small and thus we can neglect the effect of 
dispersive bands. In section UlI D1 we have seen that the 
topological pump as well as the coherent pump in the 
continuum need sufficiently long times T 0 ,Tn ^ ^ln 
avoid nonadiabatic excitations. The excitation probabil¬ 
ity to the excitonic states scale as ~ T^'^, Pq ^ ^ ^ 

(see eqns. (l20)l . (HD). Here, we show that the scaling 
for Pn also holds for the lattice in the perturbative 
regime. 

Firstly, we analyze the excitation probability Pex in 
the case of flux insertion. We start from the ground 
state |gs) with N = 1 {N = ?,) particle(s) at N^f, = 0 
= 2 ) flux quanta and insert one flux quantum in 
time to create a CF quasi-hole excitation |qh). Figure 
[5]shows the probability Pex for excitation of higher bands. 
Besides an oscillatory behavior with increasing duration 
T^, we confirm the expected scaling P^ ^ the 

perturbative regime. 

To analyze the excitation probability Pex for the co¬ 
herent CF pump, we start from the quasi-hole state |qh) 
at = 1 = 2) flux quantum. The coherent pump 

is coupled resonantly to the hole excitation for different 
bare Rabi-frequencies H. In Fig. [TUlwe plot pex for dif¬ 
ferent VL. The time To needed for a 7 r-pulse is determined 
by the maximal achievable particle number {N) (see inset 


10 ° 


10 -^ 

Pex 

10 -^ 


10 -° 


FIG. 9. (Color online) Shown is the excitation probability Pex 
to high energy states in the case of flux insertion. Starting 
with particle number N = 1 {N = 3) in the ground state |gs) 
at = 0 {N^ = 2) flux quanta, we insert one flux quantum 
and analyze the probability pex of being in an excited state 
for various flux insertion times T^. The dashed lines indicate 
the expected scaling from eqn. (I20II in section [III PI 
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FIG. 10. (Golor online) We show the excitation probability 
Pex to high energy states for different Rabi-frequencies Q in 
the coherent CF pump. Starting from the quasi-hole state 
with N,j, = 1 {N^ = 2) flux quantum, we coherently couple 
to the ground state with N = 2 {N = 3) particles for various 
Rabi-frequencies fl. The dashed lines indicate the expected 
scaling from section UlI PI In the inset, we extract the 7r-pulse 
time Tn when the maximal number of particles is in the sys¬ 
tem, for two different values of fl. 


Fig. [TUI) . For time Tn, we extract the excitation prob¬ 
ability Pex of being not in the ground state with N = 2 
{N = 3) particles. We find excellent agreement with the 
expected scaling in the perturbative regime. 


3. Dispersive bands - quasi-hole trapping 

As noted in the toy model section m the finite band¬ 
width leads to an intrinsic dispersion of the quasi-hole 
excitations. Without trapping the hole excitations, the 
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FIG. 11. (Color online) Shown is the number of particles (N) 
during the protocol in time t/T in a system with and without 
trapping potential ^h. The lattice size is 11 x 11 with an 
effective flux per plaquette a* = 0.1. We insert flux in time 
= 60/J*. The vr-pulse time Tn for ph = 0 (gh = IT*) is 
Tn = 100/J* (Tn = 130/J*). 



FIG. 12. (Color online) a) Low energy spectrum of the CF 
Hamiltonian (1241) with a* =0.1 and lattice size 11x11. States 
are labeled with number n. We identify three different sec¬ 
tors using additionally the density distribution of the states: 
Lowest Chern band (LChB) n = 1,..., 6, edge states (Edge) 
n = 7,..., 14 and higher Chern bands (HChBs) n > 15. b) 
Energy spectrum during insertion of one flux quantum </o 
with parameters as in a). 


coherent pump cannot replenish them efficiently. 

We numerically analyze the effect of dispersive bands 
on a square lattice. The results are shown in Fig. 1111 
where we compare the number of CFs (TV) {t) for 25 cy¬ 
cles of the protocol with and without quasi-hole trapping 
potential g^. Here, T = T^+Tq, is the time needed for one 
step of the protocol. The trapping potential improves the 
efficiency of the growing scheme already after three steps. 
Therefore, we include the quasi-hole trapping potential 
gii for the following discussions. 


/. Finite systems - edge decay 

Let us finally discuss the effects of a finite system, 
where edge states are present. As explained in the toy 
model section |TT1 during the protocol edge states will 
transport particles to high energy states. To reach a 
homogeneous particle density in the bulk of the LChB, 
absorbing boundaries can be implemented to prevent ex¬ 
citations to HChBs. 

In the case of open boundary conditions, we identify 
three different regimes in the single particle spectrum of 
the CFs, shown in Fig. [T^ . Due to the finite size of the 
system, edge states occur between the states of the LChB 
and those of the HChB. Crucially, the effective CF model 
on the lattice is not supposed to describe the dynamics of 
the high energy states correctly. To prepare a LN type 
state in the bulk of a finite systems, it is necessary to 
avoid high energy excitations. 

The free CF energy spectrum of a finite system with 
trapping potential g^ during adiabatic flux insertion is 
depicted in Fig. [T^ . Besides the creation of a hole exci¬ 
tation in the LChB, edge states occur connecting the low 
and high energy states. Due to the edge states, particles 
will be excited to the HChBs of the system during the 
protocol. 



FIG. 13. (Golor online) Particle number (ALChB), (ABdge), 
(Ahchb) in the three different energy sectors, defined in Fig. 
[UK, after each cycle in the protocol. The parameters used 
here are the same as in Fig. [TT]with trapping potential gn. 


This is shown in Fig. 1131 where we analyze the particle 
number (A^LChs), (A^Edge), (A^HChs) in the three different 
regions, defined in Fig. USK, after each step of the proto¬ 
col. As expected, in the first few steps the number of par¬ 
ticles in the LChB (A^lciib) increases. However, before 
the LChB is completely filled, edge states are populated 
and as a consequence the number of states in HChBs 
(A^HChe) starts to increase. Figure ITlh shows the popu¬ 
lation pn of the first CF states (labeled with integer n) 
after 25 steps of the protocol. 

To avoid the population of HChBs, we include absorb¬ 
ing boundaries as discussed in the toy model. Since edge 
states are localized at the boundaries of the system, the 
bulk properties of the system will only slightly be af¬ 
fected. However, for a properly chosen decay rate 7Edge 
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FIG. 14. (Color online) Population pn of the low energy CF 
states at the end of the protocol, a) Parameters as in Fig. [11] 
after 25 steps of the protocol and b) after 35 steps of the proto¬ 
col with absorbing boundary conditions (ysdge = 0.0025J*). 


achievable, in large systems with more than 10 particles. 

We believe, that our protocol can be used to grow other 
exotic states than the LN state, e.g. the Moore-Read 
Pfaffian [^. So far, we did not consider experimental 
realizations. However, ultracold gases as well as pho¬ 
tonic systems are promising candidates. Moreover, the 
efficiency of our scheme can be increased by introducing 
multiple pairs of topological and coherent pumps. 
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Appendix A: Flux Insertion in the IQHE 


on the boundary, edge excitations will be lost before they 
reach the HGhB during the protocol. In Fig. [T4 }d we show 
the population of the first GF states after 35 steps of the 
protocol. The HChBs are only very weakly populated 
due to the absorbing boundary conditions. The CF den¬ 
sity p* of the last cycle is shown in Fig. [TJ In terms of 
the CFs we reach an average CF filling factor of v* ~ 0.9 
of the LChB, which is close to the optimal value v* = 1. 


To address the problem of flux insertion, we briefly 
review the Landau Level (LL) problem. To this end, 
we use a basis, which is not typically used in standard 
textbooks. It will turn out, that this basis allows a simple 
understanding of the flux insertion process. 


1. Landau Level 


V. SUMMARY & OUTLOOK 


The LL Hamiltonian in symmetric gauge, A = 
B/2{-y,x,0), is {H = 1) 


In conclusion, we discussed a protocol which allows to 
grow topologically ordered states in interacting many- 
body systems. We explained all necessary ingredients 
using the SLBHM as a simple toy model. We showed 
that in the flat-band limit a combination of a topologi¬ 
cally protected Thouless pump [l^ , creating a local hole 
excitation, and a coherent pump, refilling the hole excita¬ 
tion, are sufficient. Moreover, we extended the protocol 
to the case of dispersive bands and finite systems with 
edges. 

Furthermore, we discussed the protocol in detail in 
both, the continuum case and the lattice case of frac¬ 
tional quantum Hall systems. In the continuum, we esti¬ 
mated the fidelity of the protocol depending on the par¬ 
ticle number N. To describe numerically large lattice 
systems with many particles, we introduced an effective 
model based on the GF description of the FQHE. This 
allows us to simulate large systems and include edge ef¬ 
fects. We showed, that a quasi-hole trapping potential 
can be used in the case of dispersive bands to maintain 
a high efficiency of the protocol. Moreover, absorbing 
boundaries are used in finite systems, where edge states 
would transport excitations to higher bands. We showed 
that in the case of dispersive bands and even in the pres¬ 
ence of edge states, a high GF filling factor i/* ~ 0.9 is 


TLo = - (p — Af' 

u 2M ’ 

1 


(Al) 




As angular momentum Lz is a conserved quantity, i.e. 
\HoAz\ = 0 , we can use the eigenbasis \nr,^r) of a two 
dimensional harmonic oscillator. Here, = 0,1,... cor¬ 
responds to the energy levels of the harmonic oscillator 
and G Z to those of the angular momentum. We ob¬ 
tain 


Lz\nr,(-r) = (-r\'nrT(!-r) (^-2) 

^ (2nr -I- \^r\ — ^r) lUri^r)- (A3) 

The energy spectrum of both, the harmonic oscillator and 
the LL, are shown in Fig. [T51 


2. Flux insertion 

In 1981 Laughlin (2^ explained the quantization of the 
Hall current using the argument of flux insertion. This 
idea can be used to create localized hole excitations in the 
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FIG. 15. (Color online) a) Schematic picture of the Landau 
Level construction using the two dimensional harmonic os¬ 
cillator basis \nr,lr)- b) In the flux insertion process, the 
angular momentum Ir is increased by one, while keeping the 
quantum number fixed. 


quantum Hall effect. The idea is as follows. After intro¬ 
ducing one flux quantum <j)Q adiabatically in the center, 
a circular electric field is generated. The Hall response 
thus leads to a radial outwards current creating a hole in 
the center of the system. As shown below, the quasi-hole 
excitation is quantized. 

To realize Laughlin’s argument, we include in eqn. 
(EH) the vectorpotential 



FIG. 16. (Color online) Coupling |k| from the lowest LL n). = 
0 to higher LLs rir for different angular momentum Ir- 


case, we use the basis discussed in appendix s By in¬ 
serting adiabatically one flux quantum = 4>q t/T^ in 
time ^ 1/wcj fhe angular momentum ir{t) increases 
by one, i.e. iriE) = •^r-(O) + 1. Therefore, starting from 
state |nr,£r(0)), we end in the state In^, ^r(O) +1). Here, 
we calculate perturbatively in the regime 1 

the scaling of the probability P^j, of exciting particles to 
higher LLs. 

The nonadiabatic coupling between different LLs 
Ur ^ n(, is 

94^ ~ '^dt\Tlr t ^r) 

= - idEnr,£r)/E = (Bl) 

In Fig. m we plot the coupling |k| from the lowest LL 
n(, = 0 to higher LLs Ur for different angular momentum 
£r- As expected, the coupling to higher LLs decreases. 

To estimate the scaling of Prj, with flux insertion time 
in the perturbative regime {uJcE)~^ -C 1, we consider 
a simple two level approximation with LLs nr,nr + l and 
constant coupling k. Starting from the state \nr,£r), we 
calculate in first order perturbation theory the probabil¬ 
ity for ending in state \nr + £,£r)- We obtain approxi¬ 
mately 


-A-fj, = 


Ht) ^ 

2Trr 


Defining a new angular momentum 


1 j ^ — -/iz 




(A4) 


(A5) 


we obtain the same structure as in eqn. dSl). However, 
by adiabatically increasing we change the angular 
momentum of the system. By inserting one flux quantum 
(f>o , we increase the angular momentum £r of all states by 
one, while the quantum number Ur stays fixed. This 
generates the spectral flow depicted in FIG. ITBb . 


Appendix B: Nonadiabatic Corrections in the Flux 
Insertion Process 

To calculate the probability P^ of exciting particles 
to higher LLs during flux insertion in the noninteracting 


P 4 ^ (1 - cosiuaT^)). (B2) 

[ujcEY 

We expect, that the scaling of P^ with T^f, in the inter¬ 
acting case to be the same as in the noninteracting case, 
when the cyclotron frequency oJc is replaced by the many- 
body gap Aln- This leads to eqn. (1201) . 


Appendix C: Nonadiabatic Corrections in the 
Coherent Pump Process 

In zeroth order in H/Aln ^ 1, the coherent pump 
couples the quasi-hole state |qh) with N particles to the 
LN state with A -|-1 particles. The Rabi-frequency is B 
and we choose a driving frequency Wc/2, resonant with 
the zero-interaction energy LN state |LN). Thereby, in 
the continuum no total angular momentum AL^ is trans¬ 
ferred. 
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FIG. 17. (Color online) Schematic overview over the coupling between states in different particle sectors in first order in H/Aln. 
Starting from the quasi-hole state |qh) of N particles or the LN state |LN) with A-|-1 particles, the particle sectors from N — 1 
to N + 2 are involved, fli denote the many-body Rabi-frequencies reduced by their Franck-Condon factors with respect to the 
bare Rabi-frequency fl. The many-body gaps > Aln to the excitonic states \N) are multiples of the LN gap Aln. 


In first order, we couple the quasi-hole state |qh) and 
the LN state |LN) to the excitonic states in different par¬ 
ticle sectors from — 1 to -I- 2, as illustrated schemat¬ 
ically in Fig. 1171 For simplicity, only one excitonic state 
\N) in each particle sector with many-body gap 
is shown. The effective Rabi-frequencies 17^, reduced by 
many-body Franck-Condon factors, are labeled by an in¬ 
dex z = 1,..., 4. 

For the model discussed in section nil we find, that 
the Franck-Condon factors 


L!2/n = (iV|6o.o|LN) =0 (Cl) 

fl3/n = {N-l\bo,ohh)=0 (C2) 


vanish. Furthermore, we calculate perturbatively the ex¬ 
citation probability Pq of excitonic states in first order 
in n/ApN “C 1. Starting from the quasi-hole state |qh), 
after a zr-pulse of duration Tq = nj2Q, we obtain ap¬ 
proximately the result in eqn. (1201) . There, the factor 
CTjv is defined as 


— 

<Jn — 








AN+1)2 


/^l 


.{N+2)2 


LN 


/A£ 


(C3) 


LN 


where the sum over j includes all states in each of the 
particle sectors with gap Aj'^^ > Aln- For different 
particle numbers N, we calculated the prefactor ctat in 

Fig.IZl 


[1] K. V. Klitzing, G. Dorda, and M. Pepper, 

Phys. Rev. Lett. 45, 494 (1980) 

[2] D. G. Tsui, H. L. Stormer, and A. G. Gossard, 

Phys. Rev. Lett. 48, 1559 (1982) 

[3] D. Arovas, J. R. Schrieffer, and F. Wilczek, 

Phys. Rev. Lett. 53, 722 (1984) 

[4] B. 1. Halperin, Phys. Rev. Lett. 52, 1583 (1984) 

[5] G. Moore and N. Read, 

Nuclear Physics B 360, 362 (1991) 

[6] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and 
S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008) 

[7] M. Aidelsburger, M. Atala, M. Lohse, 

J. T. Barreiro, B. Paredes, and 1. Bloch, 
Phys. Rev. Lett. Ill, 185301 (2013) 

[8] H. Miyake, G. A. Siviloglou, G. J. 
Kennedy, W. C. Burton, and W. Ketterle, 
Phys. Rev. Lett. Ill, 185302 (2013) 

[9] G. Jotzu, M. Messer, R. Desbuquois, M. Le- 

brat, T. Uehlinger, D. Greif, and T. Esslinger, 

Nature 515, 237 (2014) 

[10] M. Aidelsburger, M. Lohse, G. Schweizer, M. Atala, J. T. 
Barreiro, S. Nascimbene, N. R. Gooper, 1. Bloch, and 
N. Goldman, Nat Phys , (2014) 

[11] M. Hafezi, S. Mittal, J. Fan, A. Migdall, and M. Taylor, 


Nat Photon , (2013) 

[12] M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, 
D. Podolsky, F. Dreisow, S. Nolte, M. Segev, and A. Sza- 
meit. Nature 496, 196 (2013) 

[13] R. O. Umucalilar and 1. Carusotto, 
Phys. Rev. Lett. 108, 206809 (2012) 

[14] E. Kapit, M. Hafezi, and S. H. Simon, 
Phys. Rev. X 4, 031039 (2014) 

[15] E. Grusdt, E. Letscher, M. Hafezi, and M. Fleischhauer, 
Phys. Rev. Lett. 113, 155301 (2014) 

[16] D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 
M. den Nijs, Phys. Rev. Lett. 49, 405 (1982) 

[17] E. Berg, M. Levin, and E. Altman, 

Phys. Rev. Lett. 106, 110405 (2011) 

[18] P. Buonsante and A. Vezzani, 

Phys. Rev. A 70, 033608 (2004) 

[19] P. Buonsante, V. Penna, and A. Vezzani, 
Phys. Rev. A 70, 061603 (2004) 

[20] P. Buonsante and A. Vezzani, 

Phys. Rev. A 72, 013614 (2005) 

[21] D. Muth, A. Mering, and M. Fleischhauer, 
Phys. Rev. A 77, 043618 (2008) 

[22] R. B. Laughlin, Phys. Rev. B 23, 5632 (1981) 

[23] A. R. Kolovsky, F. Grusdt, and M. Fleischhauer, 
































14 


Phys. Rev. A 89, 033607 (2014) 

[24] R. Prange and S. Girvin, The Quantum Hall effect 
Graduate texts in contemporary physics (Springer- 
Verlag, 1987). 

[25] J. K. Jain, Composite fermions (Gambridge University 
Press, 2007). 

[26] N. K. Wilkin, J. M. F. Gunn, and R. A. Smith, 
Phys. Rev. Lett. 80, 2265 (1998) 

[27] N. K. Wilkin and J. M. F. Gunn, 
Phys. Rev. Lett. 84, 6 (2000) 


[28] N. R. Gooper and N. K. Wilkin, 
Phys. Rev. B 60, R16279 (1999) 

[29] N. Regnault and T. Jolicoeur, 

Phys. Rev. Lett. 91, 030402 (2003) 

[30] N. Regnault and T. Jolicoeur, 

Phys. Rev. B 69, 235309 (2004) 

[31] B. A. Bernevig and F. D. M. Haldane, 
Phys. Rev. Lett. 100, 246802 (2008) 

[32] S. G. Zhang, T. H. Hansson, and S. Kivelson, 
Phys. Rev. Lett. 62, 82 (1989) 

[33] N. Read, Phys. Rev. Lett. 62, 86 (1989) 


[34] J. K. Jain, Phys. Rev. Lett. 63, 199 (1989) 

[35] F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983) 

[36] Exact diagonalization is used to obtain the coefficients 

of the monomials in the Laughlin state rep¬ 

resented in the angular momentum basis. Another ap¬ 
proach would be to use Jack polynomials [3l| . 

[37] G. Brouder, G. Panati, M. Galandra, G. Mourougane, 
and N. Marzari, Phys. Rev. Lett. 98, 046402 (2007) 

[38] A. S. Sorensen, E. Dernier, and M. D. Lukin, 
Phys. Rev. Lett. 94, 086803 (2005) 

[39] M. Hafezi, A. S. Sprensen, E. Dernier, and M. D. Lukin, 
Phys. Rev. A 76, 023613 (2007) 

[40] Note, that the counting of quasi-hole excitations on a Geo 
buckyball lattice and the counting in the continuum on 
the sphere are equal and thus also the CF counting is 
correct. 

[41] G. Moller and N. R. Cooper, 
Phys. Rev. Lett. 103, 105303 (2009) 

[42] T. Liu, C. Repellin, B. A. Bernevig, and N. Regnault, 
Phys. Rev. B 87, 205136 (2013) 



